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Motivation 

Quantum  Computing 

Quantum  Control,  Control  of  Spin  Systems 
Control  of  Numerical  Algorithms 

Constructive  Controllability,  Motion  Planning  in  Robotics 
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Time-optimal  Factorization  Problem 

C-  G  compact  connected  Lie  group  with  Lie  Algebra  $ 

C  uj  :=  Qj}  finite  set  of  LA  generators  of  $ 

C  0+ :  ’’slow,  cost  expensive”  directions 
Q~:  ’’fast,  cheap”  directions 

C  Given  X  G  G,  define 


Tmin(X)  =  inf  {  ^2  | tf 

i 


X=  n 

finite 


Problem: 

•  Is  Tm in  <  oo  always?  Compute  Tmin! 

•  When  does  there  exist  a  finite ,  time-optimal  factorization? 
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Example 


Optimal  Condition  Numbers 

C  G  —  GL(n )  general  linear  group  of  invertible  matrices 

C  oo  :=  fij}  finite  set  of  LA  generators  of  g((n) 

Q+ :  ’’hyperbolic  Jacobi  rotations” 
fl- :  ’’standard  Jacobi  directions” 

L 

C  Given  X  e  G,  define  (k  denotes  the  condition  number) 


Problem: 


Imin(X)  — 


X 


n  nf } 

finite 


This  factorization  task  with  minimal  total  condition  number! 


c  Does  there  exists  factorization  with  better  condition  numbers  than 
for  XI 
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Intermezzo:  Lie  Groups  and  Lie  Algebras 


Example.  General  linear  group  of  invertible  n  x  n  matrices 


GL(n,  M) 


{X  e  Mnxn 


det  X  7^  0}. 


■  Definition.  A  matrix  Lie  group  is  any  subgroup  G  C  GL(n,  M)  that  is  also 

a  (locally  closed)  submanifold  of  MnXn. 
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Intermezzo:  Lie  Groups  and  Lie  Algebras 


Examples,  cont’d: 

(a)  The  real  orthogonal  group 


0(n)  :=  {X  e  Rnxn|  IIT  =  In} 


(b)  The  special  unitary  group 


SU (n)  :=  {X  e  Cnxn|  XX*  =  /n,  detX  =  1} 


(c)  The  Euclidean  group 


E(n ) 


R  p 
0  1 


R  G  0(ri),p  G  R 


The  first  two  examples  are  compact  groups,  while  the  third  is  not. 
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Intermezzo:  Lie  Groups  and  Lie  Algebras 


Definition.  A  vector  space  V  with  a  bilinear  operation  [ 
satisfying 


:  Vx  V-^V 


(i)  [x,y\  =  -[y,x 

(ii)  [x,[y,z]\  +  [z, 


X , 


y]\  +  [y,  [z,x]\  =  0  (Jacobi  Identity) 


is  called  a  Lie  Algebra. 
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Intermezzo:  Lie  Groups  and  Lie  Algebras 


C  Lie  algebras  are  the  tangent  spaces  of  Lie  groups. 

c  Theorem.  Let  G  C  GL(n,  R)  be  a  matrix  Lie  group.  Then  the 
tangent  space  $  :=  TjG  at  the  identity  matrix  is  a  Lie  algebra  with 
commutator  as  the  Lie  bracket: 


[X,  Y]  =  XY  —  YX. 
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Intermezzo:  Lie  Groups  and  Lie  Algebras 


Examples 

(a)  The  Lie  algebra  of  0{n)  is 


o(n)  :=  {Cl  e  Mnxn|  =  -Cl} 


(b)  The  Lie  algebra  of  SU(n )  is 


m{n)  :=  {Cl  e  Cnxn|  Cl*  =  =  0} 


(c)  The  Lie  algebra  of  E(n)  is 


e(n)  := 


Cl  v 

0  0 


=  —Cl,  v  G  Mn 
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Control  on  Lie  Groups 

G  Lie  Group  with  Lie  Algebra  g. 
Bilinear  control  system  on  G 


(S)  Xit) 


'  rn  \ 

Ad  +  ^ujWAj  X(t),  X(0 )  =  /, 
J=1  / 


where  Ad,Au  G  0. 

Reachable  Set  at  time  T  >  0 


7£(T)  =  {Xp  G  G|  and  s  <  T  :  X(s)  =  Xp} 


Reachable  Set 

K  =  Ut^(T) 
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Control  on  Lie  Groups 


Definition 

C  Accessibility:  The  reachable  set  7 Z(T)  has  an  interior  point 

C  Local  Controllability:  The  identity  /  E  1Z(T)  is  an  interior  point 

C  Controllability:  For  any  XF  E  G  there  exist  controls 
ni(-), um(-)  and  T  >  0  s.t.  the  solution  of  (E)  satisfies 

X(0)  =  I,X(T)  =  Xf. 
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Control  on  Lie  Groups 

Problem  1  (Accessibility) 

t  Definition  ( System  Lie  Algebra ) 

C  :=  smallest  Lie  subalgebra  of  g,  containing  Ai, Am,  A d 
Generators :  {[A,  B }  =  AB  —  BA) 

Ad,  Ai,  [Ad,  Ai],  [Au  Aj],  [Ad,  [Ai,  Aj]], ... 

C  Theorem.  (S)  is  accessible  if  and  only  if  the  system  Lie  algebra  is 
£  =  0. 
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Control  on  Lie  Groups 

C  Theorem  (Lian  et  al.  1994)  Suppose 

(i)  For  some  constant  controls  u\, um 

(^const)  X  =  {Ad  +  ^  UjAj)X 

3 

is  weakly  positively  Poisson  stable. 

(ii)  The  system  Lie  algebra  C  satisfies  C  =  g. 

Then  the  bilinear  control  system  is  controllable. 

Accessability  +  Poisson  Stability  =>  Controllability 
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Control  on  Lie  Groups 

Definition  (Poisson  Stability) 

Flow  of  (Dconst):  $  :  G  x  M  ->  G;  (z,£)  t— > 

4r  (Hconst)  is  Weakly  Positively  Poisson  Stable  if  for  all  z  e  G,  any 
neighborhood  B(z)  of  z  and  all  T  >  0,  there  exists  t  >  T  such  that 

$(Uz,t)nB(z)^9). 

Examples:  a  swing  (no  damping),  satellite  attitude,  ball  rolling  in  a  bowl. 
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Control  on  Lie  Groups 

i.  Theorem  (Jurdjevic-Sussmann)  Assume: 

(i)  There  exist  constant  controls  such  that  A d  +  JT  UjAj  lies  in  a 
compact  subalgebra  6  of 

(ii)  The  system  Lie  algebra  C  satisfies  C  = 

Then  the  system  (E)  is  controllable. 
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Control  on  Lie  Groups 


c  Corollary 

Let  G  be  a  compact  connected  Lie  group.  Then  (£)  is  controllable 
if  and  only  if 

C  =  0. 
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Control 


Time-Optimal  Control  on  Lie  Groups 
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Time-Optimal  Control  on  Lie  Groups 

General  Notation: 

w  Let  G  be  a  compact  Lie  Group  with  Lie  algebra  K  C  G  a 

compact  connected  Lie  subgroup  with  LA  6  .  Consider  the  bilinear 
control  system  on  G 

rri 

(S)  x  =  (Ad  +  ^  x(o)  =  i 

3  = 1 

with  Ad  G  0,  Ai, G  6. 

•  Assumption: 

•  E  is  controllable,  i.e.  g  =LA  generated  by  A&,  A i, 

•  6  =  LA  generated  by  Ai , 
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Time-Optimal  Control  on  Lie  Groups 

Given:  Initial  state  Xq  =  /,  Final  state  XF  e  G 

Problem  1 . Find  controls  'Ui(-), um(-)  s.t.  the  corresponding 
solution  X(t)  of  (£)  satisfies 

X(0)  =  X0,  X(T)  =  XF  for  some  T  >  0 

Problem  2. If  problem  1  has  at  least  one  solution,  then  find  a 
time-optimal  one,  i.e.  one  with  minimal  T  =  Topt(XF). 

Problem  1  is  always  solvable,  provided  (£)  is  controllable! 
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Time-Optimal  Control  on  Lie  Groups 

Fast  versus  slow  directions 

C  Ad  is  called  the  drift  term ,  A i, the  fast  directions 

C  Fact  1 .  If  Ad  =  0  and  (£)  controllable,  then  can  control  to  Xp  in 
arbitrarily  small  time:  Topt(Xp)  =  0,  always! 

C  Fact  2.  The  presence  of  drift  term  Ad  ^  0  is  responsible  for 

Topt  >  0. 

Idea:  Factor  out  fast  directions! 
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Time-Optimal  Control  on  Lie  Groups 

Quotient  System  and  Equivalence  Principle 

•  Consider  the  quotient  space 

G/K  :=  {Kg  \  g  e  G} 

of  left  co-sets  Kg,  K  =  exp(6)  Lie  Group  generated  by  fast 
controls. 

C  G/K  is  a  smooth  manifold 
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Time-Optimal  Control  on  Lie  Groups 

Example:  (NMR) 

C  For  the  NMR  Schrodinger  Equation  on  G  =  SU(2N) 

ujIlA  X,  X(0 )  =  I 
t  :=  LA  generated  by  Lffi,  ...,\H2n 

K  :=  exp(6)  compact,  connected  Lie  subgroup  of  SU( 2N), 
generated  by  exp(i tHj\t  G  M,  j  =  1,...,2 N. 

One  verifies  K  =  SU (2)  0  ...  0  SU (2) 

•  For  N  =  1  :  K  =  SU( 2)  =  G 

•  For  N  =  2  \  K  =  SU( 2)  0  517(2)  ~  50(4)  c  517(4) 
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Time-Optimal  Control  on  Lie  Groups 

Quotient  System  and  Equivalence  Principle 

C  The  quotient  system  of 

rri 

(S)  X=  (Ad  +  X(0)  =  I,  X(T)  =  Xf 

3  = 1 

is  the  control  system  on  G/K 

(E/K)  P  =  Adu{t)  ( Ad)P ;  P(0)  =  K,  P(T)  = 

Ad^(A^)  =  gAdg-1 ,  g  E  K.  The  control  functions  for  (E/iiT)  are 
arbitrary  functions  t  i— >  U(t)  E  K. 
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Time-Optimal  Control  on  Lie  Groups 


Quotient  System  and  Equivalence  Principle 


C  Theorem  (Equivalence  Principle). 

(£)  is  controllable  on  G  iff  (E/K)  is  controllable  on  G/K. 
Moreover,  the  optimal  times  on  G  and  G/K  coincide. 

Tgt(XF)  =  T°{k(KXf) 


Proof:  PhD  thesis  by  Khaneja 

C  The  optimal  time  T^/tK  has  an  interpretation  within 
Sub-Riemannian  Geometry. 
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Time-Optimal  Control  on  Lie  Groups 


Sub-Riemannian  Geometry 

C  Let  M  be  a  Riemannian  manifold,  E  C  TM  a  constant  dimensional 
subbundle  that  satisfies  the  Hormander  Condition 

For  any  p  £  M,  the  LA  of  the  sections  of  E  evaluated  in  p 
is  equal  to  TpM  (controllability  cond.) 

C  For  any  two  points  x,y  £  M,  the  Sub-Riemannian  distance  is 

d(x,  y )  :=  inf 


\d(t)\\dt  |  a(0)  =  x,a(  1)  =  y,d(t )  6  ^Q(t)| 


Example:  M  =  G/K,Ep  :=  span {kAdk~l  \  k  €  K}P,  P  £  M 

satisfies  the  Hormander  Cond.  (Equivalence  principle) 


C  NMR:  M  =  SJJ{2N)/ SU{2)  ®  ...  ®  SU( 2)  Sub-Riemannian  space 
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Time-Optimal  Control  on  Lie  Groups 

Sub-Riemannian  Geometry 

•  Theorem. 

T%tK(KXF)  =  d(K,KXF) 

Sub-Riemannian  distance 

C  Remark.  The  Sub-Riemannian  distance  d{x,y)  is  greater  than  or 
equal  the  Riemannian  distance  on  G/K\ 

d(x,y)  >  geodesic  distance  between  x,y 

C  There  is  one  case  where  these  distances  are  equal:  Riemannian 
symmetric  spaces. 
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Time-Optimal  Control  on  Lie  Groups 

Sub-Riemannian  Geometry 

C  Theorein.  If  G/K  is  a  Riemannian  Symmetric  Space,  then 

T0 pt(XF)  =  length  of  a  geodesic  in  G/K  that  connects  K  with  KXp 

C  Main  Advantage:  Riemannian  distances  (i.e.  lengths  of  geodesics) 
are  much  easier  to  compute  than  Sub-Riemannian  distances. 
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Time-Optimal  Control  on  Lie  Groups 

W  Theorem.  The  homogenous  space  G/K  is  a  Riemannian  symmetric 
space,  provided  (g,6)  is  a  Cartan-pair,  i.e.  $  is  semisimple  and 


0  =  e©p,  p:=tx 

satisfies 

[Mice,  [e,p]cp,  [p,p]ce 
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Time-Optimal  Control  on  Lie  Groups 


Riemannian  Symmetric  Spaces 

4*  SU(n)/SO(n)  is  a  Riemannian  Symmetric  Space 

4*  SU (4 )/ SU (2)  0  SU (2)  is  a  Riemannian  Symmetric  Space  (good! 
2-Spin  Case) 

4»  SU (8) / SU (2)  0  SU{2)  0  SU{2)  is  A/07"  a  Riemannian  Symmetric 
Space  (bad!) 
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Time-optimal  Factorization 

Let  G  be  a  connected,  compact  Lie  group  with  Lie  algebra  g. 

Let  K  C  G  be  a  connected  compact  subgroup  with  Lie  algebra  6. 
Let  A  G  g  be  a  drift  term  s.t.  (A,  1)l  =  0- 
Consider  the  discrete  control  System: 

(^d)  Xn+i  =  KnetnA LnXni  Xq  =  I  Kn,  Ln  G  iC,  tn  >  0. 


For  X  e  G  let  TQ%t(X)  : 


oo 


oo 


n—  1 


n = 1 
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Time-optimal  Factorization 

Problem: 

w  Is  (Ed)  controllable,  i.e.  does  T^pt(X)  <  oo  hold  for  all  IsG? 
C  Determine  the  “minimal"  time  ( X )  for  X  e  G. 
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Time-optimal  Factorization 


Generalized  Version  (multiple  drifts) 

C  G  compact  connected  Lie  group  with  LA  $ 

C  oo  :=  fij", finite  set  of  LA  generators  of  6 

4*  f}+ :  ’’slow,  cost  expensive”  directions 
’’fast,  cheap”  directions 

C  Given  X  £  G,  define 


Tmin(X)  =  inf{^|f+ 

i 


x=  Y[  } 

finite 
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Time-optimal  Factorization 


Problem 

C  Is  Tm in  <  oo  always?  Compute  Tmin! 

•  When  does  there  exist  a  finite ,  time-optimal  factorization? 
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Time-optimal  Factorization 


Example  1  (Euler  Angles) 


•  SO( 3),  to  =  {0+,0f}, 


1 

o 

0 

-1 

1 

o 

-1 

1 

o 

0^  := 

0 

0 

0 

1  l—\ 

a 

1 

0 

0 

1 

0 

o 

1 

o 

0 

1 

o 

C  Euler  Angles: 


x  =  e0in~  ee2nf  eo3n~ 


—7 r,  7 r 


•  We  will  show:  Euler  Angles  are  time-optimal  and 


^min  =  1 02  |  ^  [0,  7T 
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Time-optimal  Factorization 

Example  1  (Euler  Angles) 

•  SO{3),  u>  =  {0+,0+},  := 

•  Then  Euler  angles  are  i.g.  NOT  time-optimal: 

Tm in  <01+02  +  ^3  !  (Mittenhuber) 
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Time-optimal  Factorization 

Equivalence  Principle 

C  Let  G  be  a  connected,  compact  Lie  group  with  Lie  algebra 
•  Let  6  :=  (Ai, Am)L,  K  :=  expfi. 
t  Let  A  G  g  be  a  drift  term  such  that  (A,  6)^  =  0. 

c  Theorem. 

(a)  The  discrete  control  system  (E^)  on  G  is  controllable  and  thus 

T^(X)  <  oo 

(b)  For  any  X  e  G  the  minimal  times  T^pt(X)  =  Topt(X)  coincide, 
where  Topt (X )  is  the  minimal  time  for  the  control  problem 

rri 

x  =  +  x(o)  =  i,x(t)  =  x 

3  = 1 
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Time-optimal  Factorization 

Problem:  l.g.  time  optimal  factorizations  are  infinite 

Under  what  conditions  on  the  drift  term  A  are  they  finite ? 

Definition  [Haselgrove,  Nielsen,  Osborne]:  A  drift  term  A  is  called 
lazy ,  if  there  exists  e  >  0  such  that 

^oPt(etA)  <  t  for  all  t  G  (0, 6r).  (**) 

If  A  is  not  lazy,  we  call  it  fast. 
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Time-optimal  Factorization 

W  Theorein.  If  A  is  lazy,  there  are  no  finite,  time  optimal 
factorizations  for  any  element  X  G  G  —  K. 
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Time-optimal  Factorization 


Conjecture  1 :  There  exists  a  finite,  time  optimal  factorization  for 
all  X  e  G  iff  A  is  fast. 


Conjecture  2:  A  fast 


[A,  Ax]  =  0 


Remark:  Conjecture  2  implies  Conjecture  1 
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Computation  of  Optimal  Time 


Theorem  (Khaneja)  Let  (g,£)  be  a  Cartan  pair.  Let  A-1  be  the 
orthogonal  projection  of  A  onto  p  and  let  a  be  a  maximal  abelian 
subalgebra  of  p  that  contains  A^.  Then: 


4  Each  X  e  G  has  a  decomposition  of  the  form 

A  =  UTV  with  U,  V  e  K  and  £  e  exp  a. 
C  The  minimal  time  is  given  by 


Topt(A)  =  min  |t  >  0  (t  •  conv  W(A_L)^  D  exp  1(D)  7^  0| 


where  X  =  UXV  is  an  arbitrary  factorization  of  the  above  type  and 
W(A±)  denotes  the  Weyl  orbit  of  A^. 
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Computation  of  Optimal  Time 


Convex  hull  of  the  Weyl  Or¬ 
bit  of  a  ’’symmetric”  drift 
term  A 


Convex  hull  of  the  Weyl  Or¬ 
bit  of  an  arbitrary  A. 
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Computation  of  Optimal  Time 


Example  1 ,  cont’d: 


C  G  :=  SO(3)  and  g  :=  so(3), 


1 

o 

0 

-1 

i 

o 

-1 

1 

o 

0 

0 

0 

J  •  — 

1 

0 

0 

1 

0 

o 

1 

o 

0 

1 

o 

•  A:=ct^i+/?n2,  C:=(fi2) 

Euler  Angles:  X  =  e6,1^2e6,2^1e6>3^2,  ^  E 

•  Topt(A)  =  ct-1|02|, 

•  A  fast  f3  =  0. 


— 7T,  7T 
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Computation  of  minimal  time 

Example:  (NMR  cont'd) 

C  NMR-Schrodinger  equation  on  SU(4) 

4 

X  =  -2  +  *(0)=I, 

1=1 

where  Hd  :=  az  ®  az,  Hi  :=  I2  ®  <jx,  H2  :=  h  0  cry,  H3  :=  crx  0  I2, 
and  H\  ■=  ay  0 12. 

•  K  =  SU(2)  ®  SU(2). 

•  A  =  —2iriHd  and  a  :=  i(<rx  ®  ax,ay  ®  ay,az  0  <rz). 
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Computation  of  minimal  time 


Example:  (NMR  cont'd) 

Theorem.  For  all  X  =  UXV  G  SU (4)  and  U,V  G  K,  and  E  G  exp  a  fixed 
it  holds 


•  T(X) 


f  3 

=  mill  2Z  Fn 

^  n=  1 


eti27ri(crcc(g)crcc)et227ri(cry0cry)et327ri(cr2(g)cr2)  _  ^ 


} 


•  T(X)  <  | 
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Computation  of  minimal  time 

Optimization  Algorithm  (NMR  cont'd) 

Let  X(t,u)  =  U{ui,  ...,u6)Y;(ti,t2,t3)V(u7,  ...,u12), 


U(u  1  Uq)  =  0_i27rui^i0_i27ru2^2e-i27rn3i^ie-i27m4i^3e-i27m5//4e-i27m6//3 

V(uj  U12)  —  0  ^TTUj  H 1  0  i27TUgH2  0  i2'7r'Ug  H 1  0  l27TU±Q  H3  0  \  H  4. 0  \2ttu\2  H3 

_  eti27ri(crx®crx)  et22iri(cTy®(jy)  et32iri(crz®crz) 


To  compute  the  minimal  time  T(X),  we  combine  simulated  annealing 
with  gradient  methods  to  solve  the  nonlinear  optimization  problem: 

min  f(t,u)  :=\ti\  +  \t2\  +  \t3\, 
subject  to  g(t,u)  4  —  Retr(XjjX ( t,u ))  =  0 


where  t  =  [h,t2,t3],u  =  [u1:  u2, ul2]  G 


12 
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Computation  of  Time-optimal  Pulse  Sequences 


Consists  of  two  sub-problems: 

•  Given  T  >  0,  solve 

f(t,u)  <  T, 

t  >  0. 

c  Let  V(T)  be  the  global  optimal  value  of  g(t,u),  associated  with  a 
given  T  >  0. 

Minimize  T 
subject  to  V(T)  =  0, 

T  >  0. 


min 

t,u 

subject  to 
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Computation  of  Time-optimal  Pulse  Sequences 


Example 


XF  =  e 


1 

0 

0 

0 


0  0 
0  1 
1  0 
0  0 


0 

0 

0 

1 


T(Xf)  =  1.499996 

t  =  [0.499993  |  0.500017  |  0.499986] 


